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Abstract 

We study exact tunneling solutions in scalar field theory for potential barriers 
composed of linear or quadratic patches. We analytically continue our solutions to 
imaginary Euclidean radius in order to study the profile of the scalar field inside the 
growing bubble. We find that generally there is a non-trivial profile of the scalar 
field, generating a stress-energy tensor, that depending on the form of the potential, 
can be a candidate for dark energy. 
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1 Introduction 



In quantum field theory containing scalars it may occur that there are more than one 
local minima in configuration space. In such system trapped in a metastable 

vacuum will decay towards a vacuum with lower energy through quantum tunneling. 
Semiclassical methods can be used to describe the procedure [1], [2]. In this approach 
the theory is studied in Euclidean space and a classical tunneling solution that matches 
the appropriate boundary conditions is constructed. Typically this solution describes 
a bubble, which separates the true vacuum from the false vacuum, and which after its 
emission starts expanding, asymptotically with the speed of light. 

Typically it is very difficult to find analytic solutions for any given potential that 
contains a metastable vacuum. The original papers [T], [2] focus on a limit, in which 
the energy difference between the true and false vacua is small. In this limit the radius 
of the emitted bubble is large in comparison to its width, thus the name "Thin Wall 
Approximation" for this approach. This limit allows for analytical expressions of the 
bubble emission rates. However such a limit destroys all other potentially interesting 
features of the solution, especially in the interior of the bubble. 

In order to describe potential barriers that are not appropriate for the thin wall ap- 
proximation, we can either make a numerical computation, or approximate the potential 
with another one that is exactly solvable. The latter is analyzed in [3J, where a linear and 
a rectangular potential are analytically solved and the relevant decay rates calculated. 
However, if the actual potential is smooth, such potentials are clearly not a very good 
approximation for the regions of the two vacua and the top of the barrier, so several 
qualitative features of the solution may be lost. For example, the discontinuity of the 
potential in the rectangular approximation removes all dynamics as the field rolls towards 
the true vacuum. 

In this paper, we will extend the study of the triangular model and moreover solve 
some more realistic, still analytically solvable potentials, and try to extract qualitative 
features of the related physics, and possible cosmological implications. It is going to turn 
out that the triangular model with parameters of Planck scale can provide an elegant 
explanation for the order of magnitude of the measured dark energy density in our uni- 
verse. Several other options with more singular potential also exist predicting the correct 
order of magnitude. 
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2 Framework 



We are going to study exact tunneling solutions in the simple case of a single scalar 
field and a potential containing a unique false vacuum. We are particularly interested in 
analytically continuing our solutions to imaginary Euclidean radius in order to discover 
the evolution of the field inside the bubble. The Lagrangian describing our system is 

C=\{d^f-V (0). (2.1) 

From now on, we name the positions of the true and false vacua as and 0+ 
respectively and the relevant values of the potential V- and V+. 

It has been proven that in scalar field theory the spherically symmetric solutions are 
favored [I]. This is intuitively reasonable, as the sphere maximizes the ratio of volume 
to surface. Thus we assume that the tunneling solution depends only on the Euclidean 
radius p. 

Under the assumption that the solution depends only on the Euclidean radius, the 
Euclidean field equation is reduced to 

4> + -<j> = V (</>), (2.2) 
P 

where the dot implies differentiation with respect to p and the prime implies differen- 
tiation with respect to the field (p. The solution has to obey the following boundary 
conditions 

lim <j>{p) =<f> + , 0(0) = 0. (2.3) 

p— >oo 

In the following we are going to make simple assumptions for the form of the potential, 
that are going to allow us to find analytic solutions. 

2.1 Two Kinds of Solutions 

Before we proceed to solve the equation of motion, we would like to make a comment on 
the general form of the solutions. Typically we are going to assume that the potential is 
described by different formulas before and after the top of the barrier. One would expect 
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that we would result in a solution of the form 




p < R 



R- < p < R T 



(2.4) 



R T < p < R+ 



P> R+ 



This describes a bubble whom profile is as following. Outside a certain radius i?+ the 
field rests in the false vacuum. Inside this radius the field climbs the barrier between 
-R+ and Rt and then rolls down to the true vacuum between Rt and R- and then stays 
there. In such cases the analytic continuation to imaginary Euclidean radius is trivially 
(/) = </>_. In [3] we see that in the rectangular approximation, tunneling solutions always 
look like that, however in the triangular approximation we may get such a solution or 
not depending on the parameters of the potential. 

However in the triangular approximation, and as we will show later on in other cases, 
it may be true that the field never reaches the true vacuum in Euclidean space. The 
conditions leading to such a result in the triangular approximation are not very restricting 
on the parameters of the potential. In such cases the solution is going to be of the form 



and the analytic continuation to imaginary proper time is nontrivial. We will see that 
actually the field in such cases never reaches the true vacuum, but performs a damped 
oscillation around it. The rest of the paper focuses in this category of solutions. 

3 A Volcanic Potential and the Field in the Interior of the 
Bubble 

3.1 The Approximation 

The only selections of potential that preserve the linearity of the equation of motion is 
a linear and a quadratic one. The quadratic naturally is the most obvious selection to 
approximate the region of a vacuum, thus we will start our analysis studying a potential 
barrier built out of quadratics. The simplest possible barrier potential built by quadratics 



01 (P) 
4>2 (p) 

<f>+, 

V 



p < Rt 



P>R+ 



R T < p < R+ 



(2.5) 
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is described by 




Figure 1: The volcanic approximation 

and looks like in figure [TJ Because of the shape of such a potential, we call this the 
volcanic approximation. 

3.2 The Instanton Solution 

In order to find the tunneling solution we need to find the general solution to the equation 



3 I 2 
i -| — (p = m 

P 



If we make the substitution 



the equation is written as 



P 



(3.2) 



(3.3) 



p y + py - {m p + l) y = 



(3.4) 
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which is exactly the modified Bessel equation, x 2 y" + xy' — (x 2 + a 2 ) y = for x = mp 
and a = 1. Thus the general solution is 



c\I\ (mp) + c 2 Ki (mp) 



(3.5) 



where / is the modified Bessel function of the first kind and K the modified Bessel 
function of the second kind. In the construction of the instanton solution we are going 
to need the derivative of the solution. This is given by 



<p = m 



c\l 2 (mp) - c 2 K 2 (mp) 



(3.6) 



Using the above result it is clear that a solution that does not reach the true vacuum 
in Euclidean space will look like 



<t> = { 



ci-h (m-p) + c 2 -K\ (m_p) 

0- + , p<R T 

P 

c l+ h (m+p) + c 2 +Ki (m + p) 
(f>+ H , R T < p < 



(3.7) 



p>R^ 



Let's now apply the boundary and matching conditions to specify the undetermined 
constants and the radii. The solution has to be stationary at the origin. As K\ diverges 
at the origin, and 1\ is stationary, we get 



c 2 _ = 0. 



(3.8) 



Demanding continuity of the solution and its derivative at p = R + gives us the following 
two equations 



c 1+ h (m + R + ) + c 2+ Ki (m + R + ) = 0, 
c 1+ I 2 (m + R + ) - c 2+ K 2 (m + R + ) = 0. 



(3.9) 
(3.10) 



As both modified Bessel functions of the first and second kind are positive, the only 
solution to this problem for any finite R + is c\ + = c 2+ = 0. However as modified Bessel 
functions of the second kind decrease exponentially at infinity, we have the option that 
actually R + is infinite and 

ci+ = 0. (3.11) 
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Demanding that lim (ft (p) = lim (ft (p) = (ftx gives us 
p— >Rt~ p— >Rt + 

a- = - /:~^ t ^ t, (3.12) 

h (m-R T ) 

<h ~ <ft+ „ , Q 1 Q s 

C2 + = 7^~7 tt \Rti 3.13 

ifl (m + R T ) 

We have expressed all parameters in terms of Rt- Finally demanding continuity of the 
derivative at p = Rt, specifies this. 

Ki {m + R T ) h (m-R T ) = m + (ft T - (ft+ ^ ^ 

K 2 (iti+Rt) I\ {m_Rj<) m_ <^>_ — cftx 

This equation is not analytically solvable, so we cannot acquire an analytic expression 
for Rt- 

To sum up the tunneling solution is 



R T ((ft- - 4>t) h (m-p) 

T7 p < Rt 

p h (m-R T ) 

Rt (4>t - <ft+) Ki (m + p) 
p Ki {m + R T ) 



(3.15) 



where Rt is given by (|3.14[) . 

3.3 Condition for a Dumping Instanton 

We expect that in analogy to the triangular approximation |3] , if equation (|3.14|) does not 
have a solution, the solution reaches the true vacuum in Euclidean space. Both functions 
and jj^j are monotonous positive and take values between zero and one, as one 
can see in figure El 

Thus there is always exactly one solution, as long as 

m + (ft T - (ft+ , , 

a 7~ < 1 ' ( 3 - 16 ) 

m_ 0_ — <py 

otherwise we should expect a solution that reaches the true vacuum in Euclidean space. 
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Z 2 (x) 




Figure 2: The numerical part of the solution 



However, such a solution will look like 



p < R- 

ci-h {m-p) + c%-K\ (m_p) 

H , R- < p < Rt 

P 

ci+h {m + p) + c 2+ Ki (m+p) 
H , Rt < P < R+ 



(/>+, p> R+ 

Demanding continuity and smoothness at p = R-, gives us 



(3.17) 



cx-.Ii (m_i?_) + c 2 -Kx = 0, 

c X -h (m-R-) - c 2 -K 2 = 0. 



(3.18) 
(3.19) 



As the modified Bessel function are positive, the only solution to the above system of 
equations for any i2_ > is C]_ = c 2 - = 0. The only way to save this is to set R- to 
zero. Then smoothness at p = i?_ gives us c 2 - = 0. However then the field has not 
reached the true vacuum at p = R , as lim = J ^ 0. 

Thus it looks like finding a solution that reaches the true vacuum in Euclidean space 
is problematic. Indeed we can see that condition (13.160 is always satisfied. From the 
expression of the potential (|3.ip we can find 



mini- 



m. 



-){V + -V-) 



m+ 



rri- 



(3.20) 
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Thus 



77l_|_ (pT — <p— 
111- (j) + — 4>T 



m+ m l (0_ _ 0+ ) _ ^ m 2 +m 2 _ 0+) 2 _ 2 ( m 2_ _ m 2_) (y+ _ y ) 



) + Jm 2 + m 2 _ (<j>_ - (f) + Y - 2 (m\ - m?_) (V+ - V- 



, (3.21) 



which implies 

m_|_ (j)T — ' 



111- (f) + — <f>T 



m\{V + -V-)+m? + m 2 _ (<f>--(f> + )^rn 2 + m 2 _ (<f>--<f>+) 2 -2(m 2 + -m 2 _ ){V+-V-) 
" ™+ m- (m2_+m 2 _)( 9 i_- 9 i + ) 2 +2m2_m 2 _ (V+-V-) ^ 



m 2 + m 2 _(m 2 + +m 2 _){$--$ + ) 2 +2m 2 + m 2 _{V+-V-) 



That means that condition (|3.16j) always holds. Thus in the volcanic approximation, 
there is always exactly one tunneling solution, that never reaches the true vacuum in 
Euclidean space. This implies that we should expect to find some kind of damped os- 
cillation around the true vacuum in the interior of the bubble, if we analytically expand 
to imaginary Euclidean radius for any parameters of the volcanic potential. As this be- 
havior is determined by the potential at the region of the true vacuum, the fact that 
the solution in the case of the triangular or rectangular approximation may be constant 
for imaginary Euclidean distance [3], is an effect because of the non-smoothness of the 
potential. We expect that any smooth potential produces tunneling solution with the 
characteristic behavior of the damped oscillation around the true vacuum in the interior 
of the bubble. 

3.4 The Analytical Continuation to Lorentzian Spacetime 

The solution of the volcanic potential is very easy to analytically continue to imaginary 
proper time. The modified Bessel functions are analytic functions having the property 
I\ (x) = —iJ\ (ix), where J is the Bessel function of the first kind. Thus for p = it we 

get 

x , Rt (<t>- - (fr) Ji (m-T) 
9 (t) = (p- p-r, 3.23 

which clearly describes a damped oscillation of the field around the true vacuum in the 
interior of the bubble. For large r we can use the asymptotic formula for the Bessel 
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Figure 3: The Instanton Solution 



3.5 The Stress-Energy Tensor in the Interior of the Bubble 

For possible applications to cosmology, it would be interesting to calculate the stress- 
energy tensor in the interior of the bubble. Unlike traditional treatment, where the stress 
energy tensor contains only the vacuum energy of the true vacuum, here the damped 
oscillation around it is going to contribute too. 
The stress-energy tensor is given by 

dL 

T '' = em dA - at - (3 - 25) 

Applying our Lagrangian we get 

T^ v = &*<l>d v <j> - £5%. (3.26) 
Our solution depends only on proper time, thus 

d^ = Hr)^==- (3.27) 



10 



Using the above we get 

T P v = i> 2S ^L-Cfifi. (3.28) 

X\X A 

The first term of the stress-energy tensor is the kinetic energy of the remnant of the 
bubble wall inside the bubble. Far away from the bubble wall, it contributes only to the 
time-time component of the stress-energy tensor, thus it behaves as dark matter. The 
second term is identical to vacuum energy, that does not only originate from the energy 
of the true vacuum, but also on the oscillation of the field. 

Pvac = \ (-4? + ml (0 - </>_) 2 ) + V-, (3.29) 
Pwaii = 4> 2 - (3.30) 



Using our solution we get 

m 2 _R 2 T (0_ - (j) T ^ 



PvaC nr / n \ 2 

\2 



2Ji (m-R T y 
m 2 _R\ {(/>- - 4> T ) 2 



'lL(-J 1 (m„T) 2 + J 2 (m_r) 2 ) +U_, (3.31) 



Pwall = ——. 2 - J2 (m_r) . (3.32) 

L\ [vh-Rt) t 1 

For large r we can use the asymptotic formula for the Bessel function to get 



m_R 2 T (<f>_ - (j) T ) 2 

-nil (iji^Rt) 2 t 3 
2m-Rl {4>- - (t> T ) 2 



■ ( 7r \ 2 ( 7T \'' 

sin yn_r — —J + cos yn_T — —J 



+ U_, (3.33) 



P'mll - . , ' „ ,2 ■■ ' [ m - T 



j) • (3.34) 



As expected from virial theorem, the average of kinetic and potential energy terms in the 
Lagrangian cancel out, because of the quadratic form of the potential around the true 
vacuum. However this kind of ideas can in general provide a cosmological constant that 
depends on the size of the bubble. It could be a promising candidate for the explanation 
of the small size of the observed cosmological constant. 

3.6 The Decay Rate 

The volcano potential may be used to approximate decay rates of metastable vacua, in 
the case of a potential that does not fit the requirements for the thin wall approximation. 
Thus, it is interesting to calculate the decay rate. The decay rate per unit volume is 
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given by 



£ = Ae-T [1 + (ft)] , 



where coefficient B is given by the Euclidean action of the tunneling solution 



B = S E [4>{p)]-S E [<f>+] 



Spherical symmetry of the tunneling solution implies 

S E [</> (p)} = 2n 2 J 5pp 3 (-<}> (p) 2 + V[<P (p)} 
Using the formula of our solution we get 



B = S E [<t> (p)} - S E [</>+] = 2vr 2 (V- - V+) / 

•/ o 



+ 



2 r R T 



+ 



h (m— Rt) 

7r 2 n4 i$ (</>t - </>+) 2 / x 



5pp h (m^p) 2 + I 2 (m_p) 2 



Ki (iti+Rt) Jr. t 
Applying properties of modified Bessel functions we find 

vr 2 (V- - V+) 

2 h (m-R T ) . 2 



dpp K\ (m-p) 2 + K 2 (m_p) 2 



B 



2 

+ 7r 2 m_i?y 



Jl {ui-Rt) 
Finally using equation (|3.14|) we find 



„ vr 2 (V- - V+) Ri 2 o 



+ ir z m + Rr {<j>r 



K 2 {m+Rr) 
K\ (tti+Rt) 



K 2 (m + R T ) 
K x (m + R T ) ■ 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



The first term corresponds to the vacuum energy gained in the volume of the bubble, 
while the second term corresponds to the energy spent on the bubble wall. 

It would be interesting to study whether the decay rate we just calculated reduces 
to the known formula in the thin wall limit. In the limit where the energy difference 
between the true and false vacuum is small, we can derive from the form of the potential 
that 

m_</>_ + m + (f) + | e ^ f 2 \ 



+ 



m__77l_i 



-) 



+ O (e 2 ) , 



(3.41) 
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where e = V+ — V— We can use the above to find 

^±^4± = 1 £ - , + O (s 2 ) . (3.42) 

That means that the right hand side of equation (|3.14p is very close to one, thus Rt has 
to be large, as expected for the thin wall limit. This allows us to use asymptotic formulas 
for the modified Bessel functions to approximate the left hand side of (13.14p as 

* {m * RT) h {rn+RT) - 1 - -L. + O ( JjrY (3.43) 



K 2 (rn+R T ) h (m+R T ) 2/xi? T \R^ 

where 



*=-=^. (3.44) 



This allows us to calculate the size of the emitted bubble to be 



R T = 3 ^-^ +)2 + O (e°) (3.45) 



and finally substituting to formula fj3.40f) we result in 



B= W V ^ + /i|, ( , 46) 



32e 3 \e 2 
In thin wall approximation as described in pQ, the B factor equals 



B-Zg. (3.-) 



Si = / d<j)yj2V {(f)) -V-. (3.48) 



where 



In our case it is not difficult to calculate Si, when the vacua energies are close 

m_<£_+m_|_ </>_(_ 
/ m_4-mi I' 

Si = / #m+ {4> - <t> + ) + / +m+ ^ + #m_ (0_ - (fy + (e) (3.49) 



J0+ 

and after some algebra 



S 1 = M( *--* +)2 +Q( £ ). (3.50) 



Substituting the latter to (I3.47P gives us exactly the same result as (I3.46|) . 
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Another interesting limit to check is the tunneling without barrier limit, which is 
discussed in [5]. In our case this clearly corresponds to the limit m + — > 0. In this limit, 
we can use asymptotic formulas for the modified Bessel functions of the second kind to 
approximate 



K\ {tti-Rt) rri- Rt 



K 2 (m^R T ) 2 
This lets us write equation (|3,14p as 



+ 0(mi 



2 I x {m„R T ) <\>--<\>t 
The same asymptotic expansions allow us to write equation (|3,4U|) as 

b = * 2 ( v -- v +) R T + 27r2 ^2 w _ _ ^ {4>T _ M 



(3.51) 



(3.52) 



(3.53) 



Now we distinguish two cases. If <j)T — <t ) + "C 4>-—4>t we can use asymptotic expansions 
of modified bessel functions for small arguments and equation (|3.14p can be written as 



(m-ifr) 2 {m-R T f ( 6 ^ 

8 192 v ~ TJ d>_ 



(3.54) 



We can solve for Rt 



r\ 



+ 



m_ 4>- — (pT 3m_ 



<h - 0+ 



+ o 



(3.55) 



Substituting in ()3.40p and using V- — V + = \rr? 



B 



167T 2 {4> T ~ 4>+f (<P- ~ 4>T) 

3(V--V + ) 



+ o 



4>t) we find 



(3.56) 



This agrees with the results in [5] up to a factor of 2 that occurs because in our case the 
rolling region of the potential is quadratic instead of linear. 

If 4>t — 4> + ^> 4>- — 4>T we can use asymptotic expansions of modified bessel functions 
for large arguments and equation ()3.14j) can be written as 



m_RT 



OU_R%)=^^±. 



(3.57) 
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Again we solve for Rt 



R 



T ~~~2 



m_ \4> 



O 



0_ - T 



(3.58) 



and substitute in (|3.40j) to find 



B 



2tt 2 (6t - 0+) 4 



(V- - V + 
which agrees with the results of [5]. 



+ o 



b T - <j) + y 



(3.59) 



4 A Smooth Quadratic Potential and the Size of the Emit- 
ted Bubble 

4.1 The Approximation 

The volcano approximation makes a bad non-smooth description of the barrier top. We 
can use the fact that quadratic potentials are solvable, in order to improve our approxi- 
mation, and search for new qualitative properties of the solutions that originate from the 
form of the potential at the barrier top. An appropriate approximation is 

0+) 2 + y+, 4><4>i 

- <J) T ) 2 + V T , 01 < < <f)2 ' (41) 

0_) 2 + y_, 0>0 2 

which is plotted in figure [H Such an approximation provides quite a large flexibility 
in fitting an arbitrary potential. We can select the above potential in such a way that 
matches the actual potential in positions and energies of the true and false vacuum, as 
well as the top of barrier, and moreover the curvature of the potential in these positions 
or alternatively select so that it matches the positions and energies of the vacua, as well 
as one the three aforementioned curvatures and simultaneously be smooth at <p\ and 02- 

4.2 The Instanton Solution 

First we note that at the region between <p\ and 02 the equation of motion can be solved 
in exactly the same way as we did in the volcano potential, with the only deference of 



V(cP) 



-rn\ (4> 



-—m 



T ' 



\ 2 



■m 2 , (0 
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V(<p) 




v. - 



¥1 <Pt V>2 <P- 



Figure 4: The smooth quadratic approximation 



getting the Bessel functions instead of the modified Bessel function. The equation of 
motion is 

M , (4.2) 



4> + —<fi = m 2 
P 



and the general solution can be written as 

C1J1 (mp) + c 2 Yi (mp) 
= 4>o-\ . 



(4.3) 



For later use in the matching condition we also calculate the derivative of the solution 

ci J 2 (mp) + c 2 K 2 (mp) 



-m- 



(4.4) 



Similarly to the volcano potential, we expect that there are no solutions that reach 
the true vacuum in Euclidean space. We separate two cases. In the first the solution 
does not reach inside Euclidean space and in the second it does. Here we will study 
the first it is simpler and provides the information we want. Such a solution will 

look like 

cyrh (m T p) + c 2 tYi (m T p) 



P 



p< Ri 



c 1+ h (m+p) + c 2+ Kx (m+p) 
(f)+ H , Ri < p < i? H 



(4.5) 



4>+, 



p> R^ 
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Exactly as in the volcano potential case, R+ has to be infinite, and the boundary condi- 
tions for the solution imply 



C2T = 0, 
ci+ = 0. 



Demanding that lim 4> (p) = lim 4> (p) = 4>\ gives us 



CiT 
c 2+ 



n 



Ji (ititRi 
6 X - 4> + 



K x (m+fli) 

Finally smoothness at p = R\ implies 

Ki (m + R 1 ) J 2 (m T Ri) _ m + 4>i — <j) + 
Ki (m + Ri) Ji (uitRi) tut 4>t — 4>\ 

Thus the instanton solution is given by 

R\ {4>T ~ (f>i) J\ (m T p) 



+ 



p Ji (mr-Ri)' 
Ri (0i - (j> + ) Ki (m + p) 

p Ki (m + Ri) 



p < Ri 



, P > Ri 



(4.6) 
(4.7) 



(4.8) 
(4.9) 



(4.10) 



(4.11) 



where R± is given by (|4.10p . 



4.3 Uniqueness of the solution and the Radius of the Emitted Bubble 

As in previous cases the existence or non-existence of a solution is decided by the last 
equation, that occurs by the demand of smoothness of the solution, in our case equation 
(|4.10p . In figure [5] we plot and j^fy- The first graph implies that actually we are 

going to have infinite solutions. 

In figure [6] we plot the first solutions. Actually any other solution except for the first 
one is not valid, since for some region of p with p < R\ region, it is true that <f) < 4>\. 
Let' sketch a proof for that. 

Ji (m T p) 

The function describes an oscillation whose amplitude decreases monotoni- 

P. _ 

cally, as it can be seen in figure[6l Based on that and using properties of Bessel functions, 
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X 



Figure 5: The numerical part of the solution 

<p 




p 



Figure 6: The first solutions in the smooth quadratic approximation 



one can show that if ps is a stationary point of this function, then 



Ji (m T ps) 



PS 



> 



Ji (m T p) 



P 



(4.12) 



for any p > p$. Applying this for p = Ri in equation 14.111 it is easy to find that for any 
stationary point ps, with ps < R\ it is true that 



\4>s - 4>t\ > 4>t - <j>v 



(4.13) 



For any solution except for the first one there are at least two stationary points in the 
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region p < R\. For at least one of them <pg is going to be smaller than (f>T, thus making 
the above equation 

<f>S<<k' ( 4 - 14 ) 

This means that the solution for p < R\ has entered a region of values for which it is not 
valid. Thus any solution except for the first one is not valid, the solution is unique. 

This is not a general proof of the uniqueness of the solution for a general potential, 
however in this case the oscillatory behavior of the solution in the region of the barrier top 
resembles exactly the general behavior that could produce multiple solutions in general. 

So the only solution it may hold is the first solution, which is valid if </> (0) < 02, or 

else 

2Ji (m T Ri) 4>t-4>i (A 1K> 

VTLfR\ (p2 — 4>T 

If this condition does not hold we should search for a solution of the second case. 

Figure [5] implies that the curvature of the potential at the top is strongly related with 
the size of the emitted bubble. As we stated above, the only valid solution is the one 
corresponding to the smallest solution of equation (|4.10p . From figure [5] we can see, that 
this solution always satisfies 

^< Rl <^, (4.16) 
rriT rriT 

where p a ^ n is the n'th root of J a . We note that pi t i ~ 3.83 and p2,i — 5.14. So the 
characteristic radius of the emitted bubble is always of the order — . 

As this case of solutions is the one that actually lies as far as possible from the thin 
wall approximation, in which the radius of the bubble tends to infinity, we expect that 
actually serves as a general lower bound for the radius of the emitted bubble. 

4.4 The Decay Rate 

As we did in the volcanic approximation, we can calculate the B factor of decay rate 
from the Euclidean action. Using the form of our solution we find 

„3 



rtt T 

B = S E [(t> (p)} - S E [<£+] = 2vr 2 (V T -V+) Spp : 

Jo 

+ 71 m T Rl ^ T ^ f 5pp \- Ji {m T pf + J 2 {m T pf 
J\ (mTRi) Jo L 



TT 2 m\Rl (0i - c/) + ) 2 [°° r ,2^„, s2l (A , 7 s 

+ 7T~, / 6pp K i( m ~P) +K 2 (m^p) . (4.17) 

Ai (m+Ri ) Jrt L J 
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Applying properties of Bessel functions and modified Bessel functions we find 
vr 2 (V T -V + )Rf 



B 



ir m T R x {(t> T - (pi) -j— f =-r + vr m+i?! - 0+ — — . 4.18 



Finally using equation (|4,10p we find 



B 



7T 2 (V T -V+)Rf 



+ vr 2 m + i?? {<t> T - 4>+) {<Pi - 4>+) 



K 2 (m+Rx) 



(4.19) 



2 -t--xv^ ^ (m+it^ 

4.5 The Second Case of Solutions on the Smooth Quadratic Potential 

An instanton belonging in the second case will be of the form 

ci_Ji (m_/o) 



+ 



P 



p< R 2 



c 1T Ji {m T p) + c 2T Y x (m T p) 
4>t H ; , R 2 < p < Ri 



(4.20) 



P 



+ 



c 2 +Ki (m+p) 
P 



P > -Ri 



where we have already embodied the necessary boundary conditions. Demanding that 
4> (Ri) = <p\ and (ft (R 2 ) = 4> 2 and continuity results in 



ci- 



C2+ 



0i - 4>+ 



R2, 
Ri 



and 



Kxim+Rt) 

Y x (m T R 2 ) (ct> T - <fo) gi - Yi (mrgi) (fo - 0r) ^2 

Ji (ititRi) Yi (m T R 2 ) - J\ (rriTR 2 ) Yx (mr-Ri) 
■h (mr-Ri) (02 - 0r) -R2 - Ji (mTR 2 ) (4>t - 4>i) Ri 
Ji (mrRi) Yi (m T R 2 ) - J\ (m T R 2 ) Y\ (m T Ri) 

Finally demanding smoothness results in the following set of equations 

m T [c\tJ 2 (m T Ri) + c 2 tY 2 (m T Ri)} = -m_ci_/ 2 (m_i? 2 ) , 
rriT [c\tJ 2 (rriTR 2 ) + c 2 tY 2 (m T R 2 )] = m + c 2+ K 2 (m+R{) . 



(4.21) 
(4.22) 



cir 

C 2 T 



(4.23) 
(4.24) 



(4.25) 
(4.26) 
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These two equations provide a solution for R\ and i?2- Unfortunately the problem is 
very complicated, and has to be numerically solved. In figure [7] we show how a solution 
to this problem looks like. 

<p 

i i 

i i 

i i 



fi - 



<P+ 

1 I 
I I 

I L 

R 2 R] 




Figure 7: The solution for the smooth quadratic potential 



5 A Triangular Potential and a Candidate for Dark Energy 
5.1 The Approximation 

We will now study the case of a potential barrier approximated by a segment of linear 
potentials. Such a potential looks like in figure This approximation has been analyzed 
in [3]. We review this derivation and then analytically continue to negative proper time. 
In the following we use the definitions 

AV± = Vr — V±, A<p ± = ±(<p T -^ ± ), (5.1) 

A ± ,f^, CS ^. (5.2, 
In other words the potential barrier is described by 

VM = { -\ - -[ (5.3) 

4>t < <fi < 4>- 
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V(y>) 




v. - 



Figure 8: The triangular approximation 



5.2 The Instanton Solution 

In order to find the appropriate solution we need to solve the equation 

v 3 ■ . 

4> + -4> = \. 
p 

It is not difficult to show that the general solution is 

8 pr 



(5.4) 



(5.5) 



Using the above, an instanton that does not reach the true vacuum in Euclidean space 
will look like 

r a_ „ c 



-p 2 + — + K_, p<R T 
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-p 2 + — + K + , R T < p < R + ■ 



P> Ra 



(5.6) 



Let's now apply the appropriate boundary and matching conditions to determine the 
constants as well as the radii Rt and R+. The field must be stationary at the origin. 
This implies that 

c_ = 0. (5.7) 
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Demanding continuity of the derivative of the field at R = R + gives us 

Then we demand continuity at R = i?+ and R = Rt and we find 

K+ = 6 



— R+ 



A. 



3 2\2 . A- 



Oltrp O 



b T + —R^ 



Continuity of the derivative at R = Rt gives us 



:i+c)r t . 



(5.8) 

(5.9) 
(5.10) 

(5.11) 



So far we have expressed all unknowns in terms of the unknown radius Rt- In order to 
determine this final unknown we demand that (ft (Rt) = 4>t- We result in 



(VTT~c-i) 2 R T . 



(5.12) 



Now we have completely determine the tunneling solution. To sum up 



A + 1 
P 



+ 



p < R T 

R t < p < R+ 

P>R+ 



(5.13) 



where the two radii are given by (|5,lip and (|5.12p and the value of the field at the origin 
4>o is given by (15.101) . 

If we want our solution to make sense it has to be that (fto < (ft— . Otherwise the field 
has already reached the true vacuum in Euclidean space. The above condition implies 
that 

A ~ > To- (5-14) 



9 t - 9+ (VTTc - l) 

If this condition holds the field never reaches the true vacuum in Euclidean space and 
thus we expect to perform a damped oscillation around the true vacuum in the interior 
of the bubble. Otherwise the field equals exactly to (ft- inside a sphere of finite radius in 
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Euclidean space, thus the analytical continuation for imaginary Euclidean time is trivially 
4> = 4>~- Such a solution is well analyzed in [3] and we will not study it here. 




5.3 The Analytical Continuation to Lorentzian Spacetime 

We are interested in studying the solution for imaginary Euclidean radius. As we com- 
mended in section [2TTI this is interesting only if condition (|5.14p holds. In such case the 
solution is given by (|5.13p . We substitute p = ir in the formula for the solution that is 
valid for p < Rt to get 

4>o + ^r 2 . (5.15) 

This grows indefinitely as r decreases. This means that at some finite r it reaches the 
true vacuum. After that point, the solution we have is not valid anymore and we need to 
find an appropriate solution for > We approximate the potential around the true 
vacuum as 

. _ . 6 < 6- 

V{4>) = { , (5.16) 

<t> > <P- 

where obviously A~ = A_. So after the r where the solution reaches the true vacuum, 
we have to fit a solution of the form 

4> = -^r 2 + ^ + K. (5.17) 

However this solution will reach a maximum and return to the true vacuum, with some 
non-vanishing derivative, thus enforcing us to fit again a solution of the form 

^=V r2+ ^ +i ^' (5 - i8) 

and so on. So we are going to get an infinite sequence of segments describing a damped 
oscillation around the true vacuum. So we need to solve the general problem of fitting a 
solution of the form (|5.17p or (|5.18p to the boundary conditions 

0(T O ) = 0_, ^(To) = 6 , (5.19) 
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or else solve the matching conditions 



cP_ = ±^ + ^ + K, (5.20) 



A T 2c 
-o 

It is not difficult to find that the solution is 



$o = ± T T - ^3. (5.21) 



C = ± Y T ° ^ (5 - 22) 

A T „ $ T t 

K = 4>_T —Tq + —. (5.23) 

Once we found the right expression we need to find the new point where the solution 
reaches the true vacuum. Demanding that 4>(T±) = <f>-, we take 

Tl = T% T (5.24) 

0(Ti) = 6i = -6 o ^. (5.25) 

Thus we can express the solution inductively as 



A + „ c 2 „+] 

T 



—T 2 H ^2 ^ ^2n+l, T 2n+ l < T < T 2n 



A „ c 2n 

T 



(5.26) 



2 + + #2n, T 2n < r < T 2n _i 



where the constants in solution are given by 

A+ 4 6 2ra T| n A- 4 ggn-l^-i , . 

c 2n +i--— i 2n — , c 2n - — J 2n _i , O^'J 

J. , ^ + rp2 , ^2n?2„ A" 2 $ 2n _iT 2n _i ,,„„> 

^2n+i = 0- + — r 2n + — , if 2n = 0- - — T 2n _ x + , (5.28) 

and the T n 's and 6 ra 's can be calculated inductively by 

T2 _ T>2 , 4$2n-l?2n-l 2 _ ^2 4<£ 2TO T 2n 

J 2n — J 2n-1 "1 > 1 2n+l ~ 1 2n yl ) (O.Zyj 

$2n = -^n-l^ 1 , <*>2„+l = -$2n7^- (5-30) 
-£2n J 2n+1 
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Finally the initial values for T and $ can be easily calculated by equation (|5.15p (which 



R T R + 



also defines cq and Kq) 



Figure 9: The Instanton Solution 



?o = p (0- - 0o) 

The solution is plotted in figure [H 

5.4 Asymptotic Behavior of the Solution 

For large r expressions (|5.29p and (|5,30p can be approximated by 



T, 



2n+l 



2n+l 



T2n — 



2$ 



2)i. 



1 + 



2$ 



2)i 



?2n 
$2n 



^2n-l + 



2$ 



2n-l 



A+ 



$2n-l 1 



2$ 



2n-l 



2)i 



A+T^ 



2n-l 



(5.31) 
(5.32) 



(5.33) 
(5.34) 



In order to avoid the alternation of signs in we find recursive relations with step two 



2n+2 



$5 



P 2n - 2 C $ 2 n 
$2 



*2n + 2C 



2)) 



T 2 , 



(5.35) 
(5.36) 
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where c = ^ + ^r- We are interested in finding the asymptotic behavior of T 2n and 
< j > 2«- We try a solution of the form 

T 2n ~ an s , (5.37) 

6 2n ~ bnK (5.38) 

The recursive relations at leading order become 

asn s ~ l = -2cbn\ (5.39) 

Un'- 1 = — , (5.40) 



which imply that 



s-t = l, (5.41) 
-as = 2cb, (5.42) 
at = 2cb. (5.43) 



The solution is 



s = -t=\, (5.44) 

a = -Acb. (5.45) 

b in undetermined and depends on the initial values for the series, or else in the parameters 
of the potential. Thus the asymptotic behavior we are looking for is 

T 2n = -4cbV^, (5.46) 

$2n = 4=- ( 5 - 47 ) 

Using the above equations we can find asymptotic expressions for any element of the 
solution. Combining them we get 

Acb 2 

$2n = ~. (5.48) 

J-2n 
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5.5 The Stress-Energy Tensor in the Interior of the Bubble 

Now we know the asymptotic form of the solution, thus we can calculate the asymptotic 
form of the stress-energy tensor in the interior of the bubble. Equation (|5.48p implies 

1 - 2 \ A 



~ ^2- (5.49) 



Now the potential is not quadratic, which means that kinetic and potential energy do 
not average at the same value. Virial theorem implies 

{V{ct>))-V- = 2^Y (5.50) 

This means that the average value of the Lagrangian does not wash out, but asymptoti- 
cally behaves as 

-(£)~4 + ^- ( 5 - 51 ) 
If we calculate the stress-energy tensor 

T% = ^^^-L8^ (5.52) 

— X\X A 

it is now going to contain except for the kinetic energy of the wall, a cosmological constant 
like term that depends on the size of the bubble like 



A~4 + ^-> ( 5 - 53 ) 
where c depends on the specific parameters of our potential. 

5.6 A Candidate for Dark Energy 

Let's now take a wild guess. Let's suppose that big bang was the emission of a bubble, 
so we are actually living inside a bubble with the size of the universe. Let' also suppose 
that the potential of the phase transition is well fit by the triangular approximation. It 
actually suffices that the potential in the region of the true vacuum is triangular. The 
natural selection of the potential parameters is that they are of Planck scale. Finally 
let's suppose that the energy of the true vacuum is exactly zero. 

Under these assumptions, we should today observe an effective cosmological constant 
of the order 

A~J^- f (5.54) 

universe 
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which obviously depends on the age of the universe. The size of the observable universe 
today is about 10 60 in Planck units. Thus the cosmological constant we should observe 
in our scenario, would be 

A ~ 10~ 120 M^. (5.55) 

This is the right order of magnitude that we measure today [B] , [7] , [5] , [U] . 

With this model we don't solve the original cosmological constant problem [10] . 
namely why the vacuum energy of SM does not gravitate. Most probably we need more 
information on quantum gravity to resolve this. However it resolves the cosmological 
constant problems that occurred after the recent measurement of it [11]. It explains 
its order of magnitude and if we assume that matter originates from the kinetic energy 
of bubble walls, then our model also explain why the matter content and dark energy 
content of our universe are of the same order of magnitude. 

The idea of relating physical constants with cosmological quantities is not new at all. 
Paul Dirac in the 1930's observed that ratios of orders of magnitudes of cosmological 
quantities are similar to ratios of orders of magnitudes concerning the fundamental inter- 
actions. Conjecturing that this cannot be a coincidence he expressed the large number 
hypothesis [12], [13], [2], according to which fundamental constants of nature, such as 
Newton's gravitational constant, depend on the age of the universe. The recent discovery 
of another large number, namely the ratio of the theoretical and observed vacuum energy 
densities, lead to similar tries to connect the energy density of the vacuum with the age 
of the universe [15], as we do in this paper. 

Additionally the idea of the dark energy originating from a scalar is not new, too. 
Quintessence models can describe the dark energy content of the universe. In this ap- 
proach the dark energy is the effect of a slow rolling scalar field instead of one that 
performs a dumped oscillation, as in our case. The subject of quintessence is quite 
broad. A nice review is given in [16] . 

One critical objection about our model is the singular form of the potential at the 
position of the true vacuum. One should be able to find a model which predicts the 
existence of such a vacuum and moreover develop quantum field theory in the region of 
such a vacuum to show that the low energy effective description is also singular. However 
linear potential can occur by brane interactions in an orbifold like in the ekpyrotic scenario 
for the big bang [T7], [18], [19]. Definitely further study is required on this field. 
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6 The Asymptotic Damped Oscillation Inside the Bubble 



6.1 The Asymptotic Solution for Potential V = a (0 — o ) n 

It is interesting that the Stress-Energy tensor in the interior of the bubble has a direct 
dependency on the size of the bubble. It is also interesting that depending on the po- 
tential this Stress-Energy tensor may represent a substance with negative pressure, thus 
providing a candidate for dark energy. Although it is impossible to find exact solutions 
for general form of the potential, we will try to calculate the asymptotic form of this 
damped oscillation for a potential of the form V = a((p — 4>o) n - In order to simplify 
things we assume that the potential is infinite for <f> < cpQ, or else is the boundary of 
the configuration space. Thus the solution gets reflected when it reaches <pQ. This may 
represent an effective field theory, where the scalar is a moduli describing geometry of 
some brane configuration living in an orbifold. 

We are interested in the damped oscillations a field performs around the area of the 
vacuum of a potential of the form 




(6.1) 



The potential is sketched in figure [TOl 



V( V ) 




Figure 10: The form of the potential 
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Assuming a spherically symmetric solution, the equation of motion is 

$ + ~j> = -V' (</>), (6.2) 
r 

where dot represents differentiation with respect to r. If we define 

Itf + V^^E, (6.3) 



the equation of motion can be written as 



3 12 



E = <j) z . (6.4) 

r 

This form of equation gives us a better point of view about how the system dissipates. 
We cannot find an exact solution to the above, but we expect that when r is large, the 
losses are small in one period, and thus we can express the above equation as 

(£) = ~<T), (6.5) 
where T = \(f> 2 ■ Now we can use the explicit form for our potential 

y = a(0-</»o) n (6.6) 
and the virial theorem. The latter instructs us that 



Using the above, equation (|6.5p can be written 



' n+2 t V ' 



whose solution trivially is 



6n 



(E) =ct ^+2. (6.9) 



c depends on the parameters of the potential. If we assume that a characteristic mass 
scale of the potential is My then the above solution behaves as 

o n — 4 

(E) ~ (6.10) 

7-n+2 
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If we are about to explain the dark energy content of the universe as this energy stored 
in the damped oscillation of the scalar field, we should check what is the appropriate value 
for the potential mass scale My. In such case (E) has to equal the dark energy density 
measured today, and r should equal the size of the observable universe R\j. Then My 
must be 

1 n + 2 -5™- 

My~(EY^—R™-\ (6.11) 
If we use (E) = 10~ 120 Mp and Ru = 10 61 / p we get the behavior shown figure [TT1 



My 




Figure 11: The scale of the potential as function of n 
Here we would like to notice three interesting cases. 

1. n = 1. This case has already been studied in section [5j The scale of the potential 
turns out to be the Planck scale, which obviously is a natural choice. Our result 
obviously agree with the results of the previous section. 

2. n = |. In this case it turns out that My = 10~ 17 M p , thus it is about WOGeV, the 
electroweak scale. It could provide a connection to other known physics. 

3. n — > 0. In this case the size of the universe does not help us at all to solve the 
hierarchy problem between the potential scale and the cosmological scale. We just 
convert the cosmological constant problem to a classical fine tuning problem as we 
need a potential with characteristic scale the cosmological one. 
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6.2 The Stress-Energy Tensor 

The stress-energy tensor is given by 



T» u = d^d u cP-C5Z. (6.12) 



Our solution depends only on proper time, thus 



d v 4> = 4>(r)^—. (6.13) 



Using the above we get 



r ^ = ^^L-£5£. (6.14) 

— X\X A 

(I* v ) = 2^ (T) - 6? «T) -(F)). (6.15) 

In the far future or in local coordinates, this is diagonal and describes an average 
energy density and pressure 

p=(T) + (V), (6.16) 

P=(T)-(V). (6.17) 

Use of virial theorem is adequate to calculate 

w ee P - = (6.18) 

which is plotted in figure [T21 

In the previous section we distinguished three interesting cases. It turns out that the 
calculated w for those is also interesting. 

1. n = 1. In the case of a linear potential it turns out that w = — ^. This is exactly 
the boundary case between an accelerating universe and a decelerating one. 

2. n = |. In this case we find w = — |. This is not the most favored value experimen- 
tally, however it corresponds to an accelerating universe. 

3. n — > 0. This gives us w = — 1 which is the most favorable value for the dark energy 
as measured today, and indistinguishable from a cosmological constant term, as 
long as the equation of state is considered. 
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w 




Figure 12: w as function of n 
6.3 Inhomogeneities in the Dark Content 

We should not forget that actually the field performs an oscillation around the vacuum. 
That means that a fair question is whether the period of this oscillation is large enough 
in order to allow for measurable phenomena. If we neglect the dissipation term, we can 
estimate the period of the oscillation, just using conservation of E. 

r- /*>+(#)* 1 ^T(^) 1 (E\^ 

T = V2 =# = - — , n \» ' _ - . 6.19 

As the period depends on E like En~2 5 for any n smaller than 2, it turns out that as 
universe grows, the period of the oscillations is getting smaller and smaller. If we use the 
experimental values for dark energy density and the radius of the universe, we result in a 
period of oscillations at our age, that is even smaller than the Planck time, thus making 
these oscillations experimentally unmeasurable. 

Having such high frequency oscillations seems peculiar, however we study a strange 
field theory, where the potential around the vacuum is not harmonic. Further study on 
such field theories is required. 
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7 Discussion 



We analyzed tunneling solutions for some potentials that can provide us with analytic 
solutions. We learned several qualitative facts about phase transitions in field theory, as 
well as we acquired some useful tools. 

The volcano and ever better the smooth quadratic potential, can provide tools for 
calculating decay rates in several problems where the thin wall approximation pQ does 
not apply. Of course we already have the tools of rectangular and triangular approxi- 
mation [3], however, as the decay rate per unit volume depends on the Euclidean action 
exponentially, a calculation of greater accuracy may be useful. 

As the volcanic potential model teaches us, it is typical for such bubble solutions, 
that the field never reaches the true vacuum in the interior of the bubble. Instead it is 
performing a damped oscillation around it, whose amplitude is some kind of function of 
time that depends on the form of the potential around the true vacuum. This phenomenon 
may be of interest depending on the form of the potential. 

The description of the barrier between the true and false vacuum by a non-convex 
potential generates questions about the uniqueness of the tunneling solution. Naively 
such a potential generates oscillatory behavior of the field at the region of the top of 
the barrier, that could result in multiple solutions. Although we don't have a proof, the 
smooth quadratic potential example shows that this is not the case. 

An interest lesson we also learn from the smooth quadratic potential example is that 
the size of the emitted bubble strongly depends and it is of the same order of magnitude of 
the inverse of the curvature of the potential at the top of the barrier or at least the latter 
serves as a lower bound for the radius of the bubble. This has an interesting implication 
in the case of an asymptotically expanding universe. Typically phase transitions in field 
theory occur at temperatures of the order of characteristic quantities of the potential 
describing the system. Masses are such quantities. Thus according to our previous 
arguments the size of emitted bubbles is going to be of the order of inverse temperature. 
Depending on the cosmological model, the radius of the universe is also a function of 
temperature, thus resulting in an upper bound for the possible number of bubbles emitted. 
This restricts the number of bubble collisions resulting in several cosmological predictions. 
Of course this is not going to be the case in a post-inflationary universe. 

The most interesting result is that, in a triangular potential, we observe an effective 
cosmological constant in the interior of the bubble, that decreases as the bubble expands, 
in such a way that its scale asymptotically equals the geometric means of the size of the 
bubble and the scale of the potential. Assuming that the parameters of the potential are of 



35 



Planck scale, and that the true vacuum energy vanishes, we make a very good prediction 
for the cosmological constant order of magnitude, which agrees with what we measure 
today. Of course this effective cosmological constant term is not the only contribution to 
the stress-energy tensor, resulting in the ratio of pressure to energy density in the interior 
of the bubble being larger than the experimentally favored value -1. However, the model 
we use contains only one scalar and is very simplistic. The idea may be useful in the 
construction of a more realistic model. 

Later we extend our analysis to study the damped oscillation for a more general 
potential. It turns out that other interesting options also exist. If the potential is 
proportional to the square root of the field, then a potential with characteristic parameters 
of the electroweak scale predicts a dark energy content of the right order of magnitude 
and an accelerating universe. 

Of course such a model for the cosmological constant predicts a vacuum energy that 
depends on the distance from the bubble wall, thus on position inside the bubble. Current 
experiments do not rule out such a dependence, thus this is an direct experimental 
prediction of our model. 

Such a time dependent cosmological constant is also predicted by quintessence models. 
The advantage of our approach is that the value of the observed cosmological constant 
is related with the size of the universe in a more direct way. 

Moreover even in the case of a quadratic potential there is still some non-trivial form 
for the stress-energy tensor that could have interesting cosmological implications in the 
expansion of the universe. In the case of more singular potentials, a connection with 
inflation could also be interesting. Finally there is the open direction of generalizing 
to different potentials, or greater number of fields. Interesting behavior like the one we 
discovered in the triangular model, may appear in different cases. We believe that the 
most interesting future direction would be the inclusion of gravity into the problem. This 
way, we would be able to directly observe the effects of the discovered phenomena in the 
evolution of the universe. 
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